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Abstract. We investigate spectral properties of a discrete random displacement 
model, a Schrodinger operator on with potential generated by randomly 

displacing finitely supported single-site terms from the points of a sublattice of . 
In particular, we characterize the upper and lower edges of the almost sure spectrum. 
For a one-dimensional model with Bernoulli distributed displacements, we can show 
that the integrated density of states has a 1/ log^-singularity at external as well as 
internal band edges. 



While the Anderson model is by far the most studied model of a random Schrodinger 
operator and still poses open problems, recent years have seen an increased interest 
in other types of random operators such as the Poisson model, the random displace- 
ment model, or Schrodinger operators involving random magnetic fields. Some recent 
references are [TOl [TTl [21 [Sj H]. While there are good physical reasons to look at these 
models, an equally strong mathematical motivation for their investigation stems from 
challenges arising due to the lack of monotonicity properties in these models. Such 
properties have been heavily used in the theory of Anderson-type models. However, 
even in the Anderson model, non-monotonicity issues arise if one considers single-site 
terms which are not sign-definite, see, e.g. [I6l [TTl |2T] . 

Absence of monotonicity requires new ideas, which, besides posing a mathematical 
challenge, often require a better understanding of physical mechanisms, typically 
in the form of a more subtle interplay between kinetic and potential energy. In 
particular, this has become apparent in the recent works [21 |3l |15] on the continuum 
random displacement model (RDM), a random Schrodinger operator of the form 



in L^(]R'^) with single-site terms q displaced randomly by vectors uj = {uJn)n&'i from 
the sites n of the lattice U^. It was found that the lack of monotonicity can be 
widely remedied by symmetry considerations as long as one assumes corresponding 
symmetry properties of the single-site potential. This has led to key insights for the 
RDM, such as a characterization of the spectral minimum, properties of the integrated 
density of states and a Wegner estimate at low energy, ultimately leading to a proof 
of localization near the bottom of the spectrum for the continuum RDM in p5]. 

Our main goal here is to provide analogues of the results in [2J and [3] for a discrete 
version of the RDM. While we largely succeed in this attempt, in some instances 
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we fall short of carrying over results from the continuum, which is mostly due to a 
well known problem for lattice operators, the lack of unique continuation properties 
which are frequently used in j2l [3] . In particular, this is the reason why a proof of 
localization for models like the one considered here, just as for the discrete Anderson 
model with singularly distributed coupling constants, is still out of reach, see the 
remarks in Section [6l 

However, this shortcoming does not affect the ID case, where we recover all the 
results from the continuum. In fact, for d = 1 we find new phenomena, not encoun- 
tered in the continuum. In particular, for a one-dimensional RDM with Bernoulli 
distributed random displacements we find a gap in the almost sure spectrum and, in 
the case of symmetric distribution, are able to investigate the behavior of the inte- 
grated density of states (IDS) at the spectral minimum and maximum as well as at 
the edges of this gap. At all these edges the IDS has a 1/ log^-singularity and, in 
particular, is not Holder continuous. 

We mention that localization at all energies for the one-dimensional discrete RDM 
has been proven in [7], where the more general setting of random word models was 
considered. This is based on showing that the Lyapunov exponent is positive at all 
but an at most finite set of critical energies, which may give rise to quantum transport 
for wave packets with energy support close to the critical energies, while it does not 
inhibit spectral localization and also leads to dynamical localization away from the 
critical energies. 

The remaining sections of this paper are structured as follows: In Section |2] we 
introduce the discrete displacement model and state all our results. Here we also 
formulate a discrete version of the fact that "bubbles tend to the corners", a result 
originally proven in [2] for the continuum RDM which provides a central tool for our 
work. Sections |3} |4] and [5] contain all proofs, with Section [3] giving results which hold 
in arbitrary dimension and Sections |4] and [5] presenting proofs of results for the one- 
dimensional model. In particular, we discuss the one-dimensional Bernoulli RDM in 
Section [5j Section [6] contains concluding remarks including several open conjectures, 
partly based on numerical observations presented there. 

2. Model and Results 

2.1. Basics. We will construct a random potential on Z"^, > 1, by randomly placing 
single-site terms supported in a rectangular box into translates of a larger box. The 
two basic boxes are 

d d 

A := Y[[l, Mi] C Z'^ and B := 6,] C Z'^ 

i=l i=l 

where 1 < 6j < Mj and 6j, Mj G N are fixed for each i & {1, . . . ,d}. As single site 
potential we choose a function g : Z*^ — ?■ M which is supported in B. We shall always 
assume the hypothesis: 

(HI) The single-site potential q is reflection symmetric in each variable 
in the sense that 

g(xi, . . . , Xj, Xi^i, . . . , Xd) (ji^Xi, . . . , Xi—i, hi Xi -\- 1, . . . , X(i) 
for each x = (x-,)^^^ and alH G {1, . . . , d}. 
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We denote the translate of g by a G Z'^ as g^, that is qa{n) = q{n — a), n E . We 
will generally require that the support of remains in A, meaning 

d 

ae A:= J][0,M,-6,]. 

i=l 

As kinetic energy operator we choose the (negative) discrete Laplacian on 
that is, for u G l^il/), 

ihou){n) = - ^ u{m), 

\m—n\=\ 

where \k\ denotes the 1-norm of a vector k (^7/. 

2.2. Bubbles Tend to the Corners. As in [2j, a key ingredient into our investi- 
gations of the RDM will be given by a property of the single-site operators h^^^ + 
on -^^(A), where h^j^ denotes the (discrete) Neumann Laplacian on A (see Section 
below for a precise definition of the Neumann Laplacian). Define 

Eo{a) ■= mincr(/ij^ + g„), a G A. (1) 

Hypothesis (HI) implies that Eq{-) is reflection symmetric on A, i.e. 

Eo{ai, . . . , ai_i, ttj, ttj+i, . . . , ad) = Eq^gi, . . . , ai_i, Mi — hi — ai, Oj+i, . . . , a^) (2) 

for all i G {!,..., rf} and a = (ai)f=i G A. Thus, -E'o(-) is determined by a G 
riiLif'^i; Mi ~ bi], where is the least integer greater than or equal to (Mj — bi)/2. 

Theorem 2.1. Fix i G {1, . . . ,d} and aj G [rj, Mj — bj], j G {1, . . . , ci} \ {i} and 
suppose (HI) holds. If either 

(i) g 7^ is sign- definite, or 

(ii) d = 1 and Eoi^a) ^ —2 for at least one a G A, 

then Eola) is strictly decreasing as a function of ai on [ri,Mi — hi]. 

This holds in each variable, meaning, in particular, that -E'o(-) attains strict minima 
in the 2'^ corners of A. The number —2 appears in (ii) as the spectral minimum of 
Hq j^ in d = 1. See Section [g] for a discussion of the relevance of assumptions (i), (ii) 
as well as for a comparison with the corresponding result in the continuum proven in 

2.3. The Discrete Displacement Model. We now construct potentials on Z*^ by 
tiling Z"^ with translates of A and placing one copy of g into each tile. More precisely, 
if 

n := A^' 

is the set of all possible displacement configurations and u = {uJk)keZ'i C fl, then 
K; : Z*^ ^ M is defined by 

Vui ■= ^ qkM+ujk^ (3) 

where kM := {k,Mi)f^^ for k = {hYi^^ and M = (Mi^^^. For each w G fi, 

hu, := ho + (4) 
defines a bounded self-adjoint operator on £^(Z'^). 
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Figure 1. The potential V^a* corresponding to the extremal configu- 
ration oj* in d = 2 depicting the clustering of neighboring single-site 
potentials. Shaded areas represent the support of copies of the single- 
site potential q. 

The main challenge in understanding the displacement model lies in its non-monotonicity 
(in form sense) in the displacement parameters Un- For example, it is not immediately 
clear which configurations minimize (or maximize) the spectrum of h^^. It is our first 
goal to answer this question. 

The family /loj is uniformly bounded in w G fl, \\hi^\\ < 2(i-|-sup„g^ k("-)l- Therefore, 

Erriir, := luf miliaih^,) > -OO, Err^nrr 



sup maxa{h^ 



< OO. 



We are concerned with the existence of displacement configurations co which are 
spectrally minimizing in the sense that 

min a{hS} = Emin 

or spectrally maximizing, 

m?&a{h^) = Emax- 

Under suitable additional assumptions, the answer to both questions will be given 
by the configuration u* = {ouDk^z^d G Q defined by 

if ki is even 
Mj — bi if kj is odd 



for i G {1, . . . ,d}. This is the periodic configuration in which clusters of 2*^ single-site 
terms are placed into adjacent corners of their supporting tiles, see Figure [Tj 

Theorem 2.2. Suppose (HI) holds. 

(a) If either, 

(i) d > 2 and q is sign-definite, or 

(ii) d = 1 and Eo{a) ^ —2 for at least one a G A, 

then the displacement configuration u* is spectrally minimizing. 

(b) Let Eola) := mmalh^j^ — g^). // either (i) d > 2 and q is sign-definite, or (ii) 
d = 1 and Eo{a) ^ —2 for at least one a G A, then uo* is spectrally maximizing. 
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2.4. The Random Displacement Model. The displacement model is called ran- 
dom if w = (wA:)fcezd G is a collection of i.i.d. A-valued random variables. That is, 
the Uk, k G Z'^, are independent, and 

F{uk = n)= /i(M), 

for a given (fixed) distribution /j, on A, meaning > for all n G A and 

= 1. P is realized as the infinite product measure, indexed by A; G Z^, 
of the measures fi on A. By E we denote the expectation with respect to P. 

In this case, the random displacement model is ergodic with respect to shifts 
{TjUj)n = oJn+j on VL and therefore, in particular, has deterministic spectrum. Thus, 
there exists a closed set S C M such that alh^) = S almost surely. In fact, we can 
characterize S in terms of the following "periodic support theorem" . This is essen- 
tially a special case of results presented in [13], needing only slight adaptations due 
to the fact that our model is ergodic with respect to a sublattice of Z'^. 

Theorem 2.3. For the random displacement model, one may take 



s= U (5) 

where 

Cper = {cu E : UJ is periodic, fi{cOk) > Vfc G Z'^}. 

Here we call a displacement configuration u periodic if cOn+jK = for all n, j G Z'^ 
and a fixed vector K = {Ki)f^i with positive integer components Ki. 

A vector n G A is called a corner of A if nj G {0, Mj — 6j} for all i G {!,..., d}. 
Throughout this paper, we make the following assumption on the distribution /j,: 

(H2) The distribution /i satisfies fi{n) > for all corners n G A. 
Applying ([s]), an easy consequence of Theorem 2.2 is the following characterization 
of inf S and sup S : 



Corollary 2.4. //(H2) and the assumptions of Theorem 2.2 (a) and (h), respectively, 
hold, then the upper and lower edges of the almost-sure spectrum S are characterized 
in terms of u* by 

infS = mincr(/i(^*) 
sup S = maxcr(/;,^.). 

2.5. A Uniqueness Result. While u* is clearly not the unique spectral minimizer 
within all u E fl (in fact, we have mia{hi^) = inf E for almost every u G fi), it 
makes sense to ask if u* is the unique periodic minimizer. In dimension one this has 
a negative answer, but we are able to characterize all periodic minimizers. This will 
be used in the proof of our results on the integrated density of states described in the 
next subsection. 

Let d = 1 and L G N, and 5*^ denote the set of all L-periodic displacement config- 
urations, i.e. 

Sl := {to G [0, Ml - bif : Uk+i = Wfc VA; G Z}. 
For a displacement configuration u G [0, Mi — bi]^ define the numbers 

n\uj) = A;G andwfc = 0} 

n^(u) = #{A; : A; G {1, . . . ,L} and Wfc = Ml - 6i}. 
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Theorem 2.5. Assume that min(T(/i^^ + g^) 7^ —2 for some a G [0, Mi — bi]. Then 
u ^ Sl is spectrally minimizing if and only if L is even and = n^{u)) = L/2. 

If, on the other hand, min(T(/i^^ — g^) 7^ —2 for some a G [0, Mi — bi], then the 
same set of configurations characterizes all periodic spectral maximizers. 

We discuss our expectation for uniqueness results in ci > 2 in Section |6] 

2.6. The Bernoulli Displacement Model. We conclude with a more detailed in- 
vestigation of a one-dimensional special case of the displacement model, which ex- 
hibits some unexpected phenomena. Here we divide Z into neighboring pairs and, 
for each pair, randomly place a single site term into one of the two points of the 
pair. With our notation from above this corresponds to d = 1, Mi = 2, 61 = 1, 
q = X6i, where A G M \ {0} is a fixed coupling constant, f{ujk = 0) = p G (0, 1), and 
P(a;fc = 1) = 1 — p. We will refer to this as the Bernoulli displacement model (BDM) 
and denote it by h^ x, keeping track of the dependence on the coupling constant. 

2.6.1. Almost Sure Spectrum. We denote by Sa the almost-sure spectrum of h^j^x- 
Theorem 2.2 applies and thus the upper and lower edges E±{X) of Sa are given by 
maxcr(/i(^.,A) and m.ma{hi^*^x), respectively, where u* corresponds to the 4-periodic 
potential with values (0, A, A, 0) in each period. By Floquet theory 

a{h^,,x) = {E eR: \DiE)\<2}, 

where the discriminant D{E), i.e. the trace of the monodromy operator, may be 
explicitly calculated as 

D{E) = E^- 2XE^ + (A^ - A)E'^ + AXE + 2 - Al 

The observation that (y{h^^*^x) is symmetric to E = A/2 suggests to substitute E = x + 
A/2, after which D{-) becomes bi-quadratic in x, which allows to explicitly determine 
the four bands of cr(/iaj*,A)- In particular, we find that 



E±(A) = ^±^2+^ + v/4TA^. 
Moreover, the central gap (G'_(A), G'+(A)) of (y{h^*^x) is given by 



G'±(A) = ^±Y^2 + ^-v/4TP. 

A deeper fact is that ((^-(A), G+(A)) is a gap of aih^^^x) for every configuration w, 
and thus also of the almost sure spectrum Sa- 

Theorem 2.6. For every A G M \ {0} and every uj eVL, 

(G_(A),G+(A))na(VA) = 0. 

The 4-periodic operator h^*^x has two additional non-trivial gaps, located sym- 
metrically to the left and right of (G_(A), G+(A)). However, at least for |A| < 2, 
these gaps are filled in entirely by spectra from other configurations. In fact, con- 
sider a{ht^i^x) for the constant configuration := (...,1,1,.. .) giving the 2-periodic 
potential (. . . , 0, A, 0, A, . . .). As it turns out, for details see [I8], for |A| < 2 the two 
bands of cr(/i^i a) fully cover the left and right gaps of cr(/i(^*,A)- Thus, combining 



Theorems 2.3, 2.6 and Corollary 2.4, we get 
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Corollary 2.7. //|A| < 2, then 

j:, = [e4X),e4X)]\{G4X),g4X)). 

If |A| > 2, then the bands of cr{h^^i^x) cover the left and right gaps of a{h^*^\) only 
partially. In Section [6] we state a conjecture on the structure of for |A| > 2. 

2.6.2. Integrated Density of States. For the one- dimensional symmetric BDM (i.e. the 
case p = 1/2), the integrated density of states (IDS) shows surprising behavior near 
band edges. A similar result, meaning in particular that the IDS is not Holder con- 
tinuous at certain energies, was first shown in an analogous setting for the continuum 
displacement model at the bottom of the spectrum in [3J. For the discrete case con- 
sidered here we get that the same phenomenon appears not only at the lower and 
upper edges of the almost sure spectrum, but also at the edges of the central gap 
identified above. 

To define the IDS, let L G N. For /c G Z, we set = [2k - 1, 2k] and define 
A(L) = uf^iAi = [1,2L]. Let hl^ ^ denote a restriction of h^^ x to A(L) with arbitrary 
boundary condition (i.e. choice of the diagonal matrix elements at —L and L). Set 

eM^,) < eM^,) < . . . , 

the eigenvalues of h^ ^ counted with multiplicity. The IDS of /i^j.a at i^^ G M is 
NxiE) = \im ^J^E(#{A; G N : E^^^) < ^}), 

which exists due to ergodicity of our model, e.g. [12], and is independent of the 
boundary condition. 

First, we note the following symmetry property which simplifies matters. 

Theorem 2.8. If E±{X) denote the upper and lower edges of the almost-sure spectrum 
of the symmetric BDM (i.e. p = 1/2) and Nx denotes the IDS, then 

Nx{E4X) +t) = l- Nx{E4X) - t) (6) 

for every t G M. 

We want to describe the asymptotics of the IDS near the four band edges -^^^(A) 
and G±{X). Due to symmetry we only need to consider the lower band edges -E-(A) 
and Gj^{X). 

Theorem 2.9. Fix A > and let Eq = E_{X) or Eq = 6*+ (A). Then there exist 
constants G,e > such that 

Nx{E) - NxiEo) > forallEe{Eo,Eo + e). 

log {E - Eq) 

A corresponding result, for energies E to the left of Eq, holds at the upper edges E^{X) 
and G_(A). 

For more discussion, including a conjecture on the asymptotics of the IDS at pos- 
sible additional band edges for the case |A| > 2, see Section [61 
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3. Bubbles Tend to the Corners and Consequences 



Our first goal in this section is to prove Theorem 2.1, i.e. that "bubbles tend to 



the corners" . This will be done in Section 13.21 after Section 13.11 will introduce the 
discrete Neumann Laplacian and list its relevant properties. The characterization of 
the spectral minimum of the displacement model, i.e. Theorem 2^ and Corollary |2.4 



3.3 



are consequences of Theorem |2.1| and will be proven in Section 

3.1. The Neumann Laplacian and Basic Properties. Let A C Z'*. The trunca- 
tion operator hQ \ on A is the operator on £^(A) with matrix elements 

-1, e A, p - j||i = 1 

0, otherwise 

The edge counting function on A is the function r^A : A — )■ N U {0} given by 

nA(0 = #{jeZ'\A:(z,j)e5A}. 

Here dA is the boundary of A, i.e. 

OA = G Z'^' X : z G A, J ^ A, IK - j||i = 1 or 2 ^ A, j G A, ||2 - = 1}. 

Associated to the edge counting function is the edge counting operator A^a specified 
by the matrix elements 

AT (, ,^ _ / M^^, i = j, ^ e A 

iU[i,3) - I otherwise ' ^' ' 

Definition 3.1. //Ac Z*^, then the d- dimensional discrete Neumann Laplacian on 
A is the operator /i^a defined on £^(A) and given by the operator sum 

The matrix elements of h^^ are given by summing the corresponding matrix elements 

of /io,A dlT'd N\. 

We now summarize some basic properties of the Neumann Laplacian. The above 
definition and most of these properties (in fact, probably all) can be found in various 
references, e.g. [201 ttH [12] ■ Detailed proofs of all of them can be found in [IB]. 

The first property is a "reflection principle" for the Neumann Laplacian, meaning 
that a solution u to the equation (/i^a + l)'^ = -^"^ can, by reflection, be used to 
construct a solution on a larger set A'. More precisely, let A = HiLif'^i^^i] 
g, M : A — )■ M with 

«A + = Eu (8) 

for some G M. Let k G {1, . . . be fixed and Uk,Tef (respectively, qk,ref) be the 
extension of u (respectively, q) to the set 

fc-l d 

A' = Y[[ai, hi] X [ak, 2bk - Ofc + 1] x JJ [a^, bi] 

i=l i=k+l 

obtained by reflecting u (respectively, q) about bk + ^ iia the kth component. It is 
not too hard to show that (/io,A' + qk,ref)uk,ref = Euk,ref- We summarize this result, 
along with three others, in 
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Proposition 3.2 (Properties of the Neumann Laplacian). 

(i) (Reflection Property) If u satisfies then 

(^0,A' + Qk,ref)Uk,ref = EUk,ref- 

(ii) (Neumann Splitting Formula) //Ai C A C Z'^, then 

'io.A ^ '''0,Ai '''0,A\Ai 5 

in the quadratic form sense. 

(ill) ( Simplicity and Positivity of Ground State ) If A C Z'^ is connected 
and g : A — )■ M, then the ground state eigenvalue of h^ j^ + q is simple and the 
corresponding eigenfunction may be taken strictly positive. 

(iv) (Ground State Energy of h^^) If A C , then mmaihJ^ = —2d. 

These properties will be used frequently in the proofs of our results. 



3.2. Proof of Theorem 2.1 , Theorem 2.1 is a statement about Eo{ai, . . . , aa) as a 



function of the ith component Oj with all other components held fixed. Without loss 
of generality, we may assume i = 1. Then since a2, . . . , are held fixed, to simplify 
notation, we write 

Qai '■= Q{ai,...,ad) 

H{ai) := hQ^^ + qa,. 
For ai G [ri. Mi — 6i]nZ, let Ua^ denote the (up to a constant multiple) unique positive 



ground state corresponding to -^0(01, 02, ... , ad) known to exist by Proposition 3.2 

H{ai)ua^ = Eo{ai, 02, ... , ad)uai. (9) 

Let M*^ and g^^ denote the extensions of and qa^ to 

A* := [l,2Mi] X [1,M2] X ■■■ X [l,Md] 

which are symmetric about the axis Xi = Mi + \. Let Ua^^ext (respectively, qa^^ext) 
denote the extension of m*^ (respectively, q^J to the strip 

Z X [1,M2] X ■■■[l,Mrf] 

which is 2Mi-periodic in the first component. 

That satisfies the ground state equation ^ implies 

Ka' + Ia*)< = ^o(ai, • • • , a,)<. (10) 

Now we turn to the proof of the theorem. Our goal is to show 

Eo(ai + 1, 02, ... , ad) < Eo{ai, 02, ... , ad) 

for all Oi G [ri. Mi — 61 — 1]. Let ai G [ri. Mi — 61 — 1], note that the former implies 
that qa-^ (1, X2, . . . , Xd) = for all a;j G [1, Mj] and i E {2, . . . ,d} which guarantees that 

= (1, X2, . . . , X,) = q:, (2Mi, X2, . . . , Xd). (11) 

If R is the right shift in the first coordinate, {Rip){ni, n2, ■ ■ ■ , rid) = "^(^i— 1, ■"-2, • • • , "^d), 
and 

fa\,R • Rfa\,ext\K*i 
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for / G {u, q}, then 

(x2,...,xa)&A 



(12) 



where 



A:= [1,M2] X ... X [l,Mrf 



We encourage the reader to check the above calculation of ( 12 ) (at least in the simplest 
case d = 1) using the definition of the Neumann Laplacian, reflection symmetry of 



g^jA* in the flrst coordinate, and (11). 



Obviously, the quantity ( [T2| ) is non-positive. If (12) vanishes, then the restriction 
of Ua^ to the "slab" K := {1} x [1,M2] x . . . x [l,Mrf] turns out to be the ground 



state of the Neumann Laplacian on K, for details see [18]. Thus, by Proposition 3.2 
Eo{ai, . . . , ad) = —2d. In other words, if Eo{ai, 02, ... , ad) 7^ —2d, then the quantity 



in (12) is strictly negative. 



The assumptions of Theorem 2.1, in either case (i) or (ii), are enough to guarantee 
that EQ{ai,a2, ■ ■ ■ , ad) 7^ —2d. In case (i), this can be seen by noting that the ground 
state eigenvalue of h^^ is —2d, and that a sign-deflnite perturbation strictly increases 
(g > 0) or decreases (g < 0) the ground state eigenvalue. In case (ii), we make use of 
the following fact: 



If min(T(/io^ + q, 



-2 for some a G [0, Mi — bi], then mincr(/iQ^ + Qa 



for all a G [0, Mi — bi]. To see this, note that mma{hQj^ + g^) = —2 implies that 
the ground state is constant outside the support of g^. Thus the ground state can be 
shifted together with the potential to produce the ground state for other values of a, 
leaving the ground state energy unaffected. 

Thus, under the assumptions of Theorem 
meaning that 



2.1 



the quantity (12) is strictly negative. 



{iK,A* +QauR)UauR,UauR) < ((/io,A* + I A* Xi , <i ) = £^0(^1, ^2, . . . , Od) |Ki 



* l|2 



where the equality follows from (|10|). 
By deflnition, |KJ| = \\ua^,R 



(13) 



splitting formula from Proposition 3.2 along with (13) gives 



applying the variational principle and using the 
37 

Eo{ai,a2, . . . ,ad) > inf a {h^^^, + qa^^n) (14) 
> mia{H{ai-l)®H{ai + l)) 

= mm{Eo{ai-l,a2,...,ad),Eo{ai + l,a2,...,ad)}. (15) 



We prove the theorem by induction on ai using (14) and (15). The flrst step is to 
show 

Eo{ri, 02, ... , ad) > Eo{ri + 1, 02, ... , ad). 
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We make use of (15) above with ai 
symmetry of Eq{ 



Ti. If Ml — hi is even then by reflection 



in the first coordinate on A, 
Eoiri - 1, 02, ... , ad) = Eo{ri + 1, a2, . . . , a^), 
thus the minimum in (15) is certainly Eo{ri + l,a2, . . . , ad). If Mi — bi is odd, then 

Eo{ri - 1, aa, . . . , a^) = Eo{ri, 02, ... , ad), 



and because we have the strict inequality (14), the minimum in (15) is -E'o('"i + 
1, 02, ... , ad). 

For the induction step, suppose Oi — 1, Oi, ai + 1 G [ri. Mi — bi] and 
Eo{ai - 1, 02, ... , ad) > Eo{ai, 02, 



,ad) 



Using (15), we have 



Eo{ai, 02, ... , ad) > min{£'o(ai - 1, 02, ... , ad), E^iai + 1, 02, ... , ad)}. 



With the induction hypothesis ( 16 ), the minimum in ( 17) must be EQ(ai+l, 02, 



(16) 

(17) 
,ad). 



3.3. Proof of Theorem 2.2 and Corollary 2.4 In addition to Theorem 2.1 the 
proof of Theorem 2^ relies on the following two facts, which are discrete versions 
of results known as AUegretto-Piepenbrink Theorem and Shnol's Theorem in the 
continuum, e.g. p!9] . 

Theorem 3.3 (AUegretto-Piepenbrink). Let E eR and V : Z'^ ^ R be bounded. If 
there exists a non-trivial u > satisfying the difference equation {Hq + V)u = Eu, 
then E < inf a{ho + V). 

Theorem 3.4 (Schnol). If u is a polynomially bounded generalized eigenfunction of 
Hq + V corresponding to E gR, then E G ct^Hq + V). 

Proofs of these facts, most likely also well known, are provided in |JJJ. 



Proof of Theorem 2^. (a) Note that restricted to A+kM with Neumann boundary 
condition is unitarily equivalent (via translation) to h^^ + g^^,. Thus, the splitting 
formula for the Neumann Laplacian from Proposition |3 2 gives 



inf a{h^ 



> 



> 



inf a( ^{h^q^,: 



inf min a{h 



is a lower bound on the 



> min {-E'o(a) : a G A} 

= EQ^amin), 

with the corner position Omm = {Mi — 6i)f=i. Thus, Eo{amin 
infimum of the spectrum of any operator /i^ with w G fi. 

In view of (18), co* is spectrally minimizing if inf a{hi^*) = EQ^amin)- Applying the 
AUegretto-Piepenbrink and Schnol Theorems, it is enough to show the existence of a 
strictly positive bounded function, ip, which satisfies h^*ip = i?o(amm)V' in the sense 
of a solution of a finite difference equation. 

so that 



Let u denote the strictly positive ground state of /i^A + 1a„ 



+ qa 



Eo{amin)u. 



(19) 
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Let U and Q denote the extensions of u and . to 



A2 = n[i'2M,: 



1=1 



which are reflection symmetric on A2 in each of the d coordinates. Looking at (19), 
we immediately have 

+Q)^ = ^o(a™„)f/. (20) 

If ip is the extension of U to all of Z*^ which is 2Mj-periodic in the ith coordinate, 
then reflection symmetry of U on A2 and (20) give 

h^^ip = Eo{amin)'4'- 

(b) To show u* is spectrally maximizing for the displacement model defined with 
single-site potential q, it suffices to show u* is spectrally minimizing for the model 
with single-site potential —q. In fact, applying spectral mapping and the unitary 
involution 

(f/0)(n) = (-l)^ti I"'l0(n), n G G f{Z'^), (21) 
we have UhoU = —Hq, and thus, for any u E Q, 

maxcr(/io + K;) = - mino-(-(/io + K;)) 

= -mma{U{-ho-V^)U) 

= -mma{ho-V^). (22) 



If CO* is spectrally minimizing for —q, then (22) gives 

max(T(/io + K;) < — min (t(/io — K)*) 
= max(T(-/io + K;*) 
= maxa{ho + Vuj*), 

for any u E Q. Therefore, u* is spectrally maximizing for q. 

That u* is spectrally minimizing for —q follows immediately from (a): If d > 2, 
sign-definiteness of q implies sign-definiteness of —q. By (a), sign-definiteness of 
—q guarantees that u* is spectrally minimizing for —q. li d = 1, the assumption 
mma{hQ^ — qa) 7^ —2 together with (a) implies u* is spectrally minimizing for —q. 

□ 

Proof of Corollary 2.4 Under the assumptions of Theorem 2.2[ a) we have 

mincr(/i^.) < mincr(/i^) 
for every u E Q. Thus mincr(if^*) < inf S. Since u* G Cper, we have by The- 



orem 2.3 that cr{hi^*) C S and, in particular, mina{hi^*) > infS. We conclude 
inf E = min a{hu]*). 

Similar, under the assumptions of Theorem 2.2[ b), it follows that 

supS = msixa{hi^*). 

□ 
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4. NON-UNIQUENESS OF THE OnE-DIMENSIONAL MINIMIZER 

4.1. The Discrete Periodic Laplacian. Suppose 5* C Z'' has the form 

d 

1=1 

for ttj, bi G Z, with < bi for all 1 < i < d. For m : 5 — ?■ C, define the function 
Mr : Z"' ^ C by 

Here L = {Li)f^i with Lj := 6j — Oj + l and fcL := {kiLi)f^^. Thus, ut is the L-periodic 
extension of u to all of Z'^. One may view ut as m defined on a (i-dimensional torus. 

The discrete Laplacian on S with periodic boundary condition, or the periodic 
Laplacian, denoted Hq ^, has domain i'^{S) and acts on a function u G £^(S') by 

{h^ gu){n) = {hoUT){n) for all n E S. 

Thus, if Ut is viewed as m on a (i- dimensional torus, then h^g is simply the discrete 
Laplacian on the torus. We summarize some properties of the periodic Laplacian for 
which proofs can be found in |18j : 

Proposition 4.1 (Properties of the Periodic Laplacian). 

(i) //V : S* — 7- M, then the ground state of h^g + V is simple and the corresponding 
eigenspace is spanned by a strictly positive eigenf unction. 

(ii) h^g > h^g in the form sense. 

(iii) If V : Z"^ ^ M. is L-periodic, then inf a{hQ + V) = inf <7{hQ^ + V\c), where C 
is any period cell ofV. 

(iv) Let u and V be functions from [a, b] (as an interval in Z) into M and E eM.. 
Ifu{a) = u{b), then {h^^^ ^ + V)u = Eu if and only if {h^^^ ^ + V)u = Eu. In 

particular, if V is reflection symmetric, then the ground states of h^^^ and 

^o[a b] coincide. 

4.2. A Preliminary Result. We begin with a result that provides a necessary cri- 
terion for a periodic configuration u to be spectrally minimizing for h^. This result 
holds for arbitrary dimension d > 1. 

Suppose G is an L-periodic configuration so that 

^i+{niLi,n2L2,...,ndLd) = ^i+riL = 

for every i G Z*^ and every n = (ni, n2, . . . , n^) G Z^. It follows that the potential V^j 
given by ^ is L-periodic, where L = LM = {LiMi)f^^ and a period cell for K; is 

d 

A = J][1,L,M,]. (24) 

i=l 

The cell A can be written as a disjoint union of translates of the basic cell A = 
nti[l, Mi], by setting /C = nti[l. ^ and defining A^ = A + (A; - 1)M for k E IC. 
Here A; — 1 is the vector with ith component ki — 1. 

Before proceeding with our preliminary result, we fix the following notations. For 
• G {N,P}, let /i* 5 denote the Neumann (• = A^) or periodic (• = P) restriction 
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of to the set S d 1/ and let EQ{h*^ g) denote the corresponding ground state 
eigenvalue. For uji G A, let EQ^uji) denote the ground state eigenvalue of h^j^ + q^^-. 
Lastly, set as before Emin = inf^^gQ min (T(/it^). 

The following result shows that in a spectrally minimizing periodic configuration 
u, every is a corner of A, i.e. it lies in 

C := {(ai, . . . ,ad) e A : ai e {0,Mi - bi} for alH = 1, . . . , d}. 

Lemma 4.2. Ifu E Q is a periodic configuration with mincr(/i(^) = Emin, then Uk E C 
for every k G Z'^. Moreover, in this case, Emin = -^o(^^a) — ^oih^^j^, where A is 
the period cell for given by (24) ■ If 'ipuj is the positive ground state for ^ with 



ground state eigenvalue Emin, then, for every k E IC, V't^Ufc is (up to normalization) 
the positive ground state for with ground state eigenvalue Eoluk) = Emin- 



Proof. By part (iii) of Proposition 4.1, Emin = in.m(7(h^) = miner (/i^^^). Part (ii) of 
the same proposition along with the variational principle gives Eo{h^ j^) > E^ih^j.^). 



Moreover, since tp^j minimizes the the quadratic form for 



> 



> E^oK)|^ (25) 

_ / , -L-'min / I IS J-'min V'^^/ 

where ip^ := ipu]\kk- The second to last inequality follows from the variational principle 
and the fact that i?o(wfc) is by definition the lowest eigenvalue of /ifJ^A^- The last 
inequalit y fol lows from the fact that Eo^Uk) > Emin as E^^Uk) > -Eo(c) = -Emm by 
Theorem 2.1 for any c E C. Division by the quantity [tpj^, ■?/'^) is permitted since strict 



positivity of implies strict positivity of ip^ and thus, {il'^,i(>^) > 0. 

We conclude that for a periodic minimizing configuration u all inequalities in the 
above chain are actually identities. 



If at least one Uk does not belong to C, then the last inequality, (26), is strict in 



light of Theorem 2A^ Thus, cok E C for all k E JC. The fact that 



> Eo{uJk) for allkEJC (27) 
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was used to obtain (25), and if at least one of the inequalities (27) is strict, then (25) 
is strict. Thus, there must be equality in (27) for every A; G /C. Therefore, ip^ is the 
ground state for /i^A^ Eo^Uk) = Emm is the corresponding lowest eigenvalue. □ 

4.3. Characterization of One-dimensional Spectrally Minimizing Configu- 
rations. Now, we specialize to the one-dimensional setting. The proof of Theorem 



2.5 relies on the following basic fact about the "shape" of a one-dimensional ground 



state eigenfunction. 

Lemma 4.3. Suppose a G [0, Mi — 6i] fl Z and a ^ ^^-^ ■ If is the strictly positive 
ground state corresponding to h^^^ M^]~^Qa with ground state eigenvalue E and E ^ —2, 
then ^ ip{M). 

Proof. By symmetry of the potential q in and uniqueness of the ground state, 

it is enough to consider a G [r. Mi — where r is the least integer strictly greater 
than (Ml — bi)/2. We may take ip to be normalized. 

Let qa,ref dcuotc the reflection of qa at the center of [1, Mi], i.e. 

qa,ref{n) =qa{Mi + l-n), r2G[l,Mi]. 

By symmetry of g on [1, bi], the normalized strictly positive ground state eigenfunction 
of Hq^^ ^^^] +qa,ref is ipref, the reflection of ip on [1, Mi], and the corresponding ground 
state eigenvalue is E. 

If = ipiMi), then it follows from part (iv) of Proposition 4.1 that the ground 
state eigenvalue of h^^^ + qa is E and ip is the corresponding eigenfunction: 

«[i,M.] + Qa)i^ = Eij. (28) 

Let Q and \E' denote the Mi-periodic extensions of qa and ip, respectively, to all of Z. 
Using the definition of the periodic Laplacian and ( 28 ) gives 

{ho + Q)^ = E^ 



which implies by Theorems 3.3 and 3.4 that Emin = inf (^{ho + Q), since \E' is strictly 
positive and bounded (thus, polynomially bounded). Therefore, by part (iii) of Propo- 
sition 4J^, E is the lowest eigenvalue of the periodic restriction of Hq + Q to any period 
cell, and \E' restricted to the same period cell is the corresponding ground state eigen- 
function. 

The trick is to choose a period cell C for which Q\c is, up to translation, identical 
to qa,ref- The appropriate period cell is 

C = [6i + 2a + l-Mi,6i + 2a]. 

Note that the assumption a > r implies that Mi -|- 1 G C. We have 

Kc + Q\cmc = E^\c 
^lc(Mi) = V'(Mi) = ^^J{l) = ^\c{Ml + l). (29) 

Since C is a period cell for \E', we also have ||\E'|c|| = Wi'W = Wi^refW- 

By construction, \E'|c must be a constant multiple of the translate of tpref to C. 
Normalization requires this constant multiple to be one, thus 

^\c{n) = ilJref{n- {bi + 2a- Ml)), n e C. 
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Therefore, we have that 

*|c(Mi + 1) 

and it fohows that 



Iprefi'^Mi - hi- 2a) 

iprefi^Mi - 6i - 2a H 



1), 



ip{bi + 2a - Ml + 1) = ij{bi + 2a - Mi] 



(30) 



(31) 



Note that 6i + 2a — Mi < 1 + a and 1 + a is the minimum of the support of qa- 



Therefore, (31) is a contradiction to the basic fact that the ground state ip is either 



strictly increasing (if E < —2) or strictly decreasing (if E > —2) below the support 
of Qa, |IHl Lemma 3.9]. 

□ 



Now we have everything we need for our proof of Theorem |2.5 



Proof of Theorem 2.5. It follows from (22) that u is spectrally maximizing for /io + K; 
if and only if it is spectrally minimizing for -V^^. Thus the claim made in The- 
2.5| for spectral maximizers follows from the result on minimizers. It therefore 



orem 



remains to prove the latter. 

Let u be an L-periodic minimizing configuration, so that En 



Lemma 4.2, it must be that uji G {0, Mi — hi\ for alH G Z and that 

E ■ 



min(T(/i^ 



miner (/i,^ 



mino"(/itj). By 



(32) 



where A is the period cell [1, LMi\. Therefore, u E Sl and n^{uj) + n^(w) = L. 

Let denote the strictly positive ground state of the periodic operator hj^^^ corre- 
sponding to Emin- Applying part (ii) of Proposition 4.1, along with the variational 
principle and (32), gives 



E„ 



(^,V^) 



> inf 

0^0 



E„ 



Thus, both inequalities above are actually equalities and ip is the ground state of 
^fj^A corresponding to Emin- From the ground state equations h^j j^ip = Emin4' for 
• G {N, P}, one gets 

-tpiLMi] 
The two together imply 



-^(2) + K.(i)V^(l) 
-^(2) + K.(i)^(l) 



E„ 

Er, 



^V^(1), 
^^(l)• 



Set Afc 



state of Hq j^ + QMi-b 



^(1) = tP{LMi) 
A + {k- l)Mi with k e {0,...,L}. 



normalized so that 



(33) 

1 denote the positive ground 
1 denote the positive 



Let 

1, and let 

ground state of h^j^ + q normalized so that 0-i(l) = 1. 

Since the potential q is reflection symmetric, given by 0(1 + n) = 0i(Mi — n) for 
n G [0, Ml — 1] is also a positive ground state for h^^ + q. Uniqueness of the positive 



ground state provides that 

1 = 



;> = 0i(Mi)0_i. Thus, 

H(l) = 0(Mi)=0i(Mi)0_i(Mi) 
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Therefore, if aii^uj) = 1 (respectively, ai{u) = — 1) w hen Ui = Mi — bi (respectively, 
Ui = 0), then (t>a,{u^){Mi) = 0i(Mi)°»('^). Lemma |43| guarantees that 0i(Mi) 7^ 1. For 
simplicity of notation, set p := 0i(Mi). 

Using the reflection property of the Neumann Laplacian, the function \E't^, con- 
structed by concatenating rescaled copies of and defined piecewise on each Afc, 
ke{l,...,L} by 

= p^'<^"'("Va.(a.)(^ - A;Mi), neAk 

(with the convention that an empty sum is zero) satisfies j^^ui = Emin^ui- There- 
fore, we have recovered, up to a constant multiple, the ground state ip. Moreover, 



P 



This is seen as follows: 



^,(A;Mi) =p^><'="'(-Va.(a.)(Mi 



P 



Since coincides at the endpoints of [l,LMi], (33), it must be that 



^,.(11 



P 



P 



nO(aj)-ni(tj) 



Since p 7^ 1, we conclude n^{uj) = ri^^u). Using n^^u) + n^{uj) = L gives 2n^{uj) = L, 



so L is even and n (u) 



L/2. 



□ 



5. The Bernoulli Displacement Model 



In this section we carry out the remaining proofs of Theorems 2^, 2^ and 2.9 
which establish properties of the Bernoulli displacement model. 

5.1. Almost Sure Spectrum. 



Proof of Theorem 2^. By spectral mapping, it suffices to show that 

i?,..A := ( /i,. A - ^ 1 - (2+ — ] > -V4 + A2 for all co e Q. 



2 J V" ' ^ 
Fix u E Q. The operator H^^ x is a five diagonal operator 



H, 



s.(-2) sU-l) 1 

1 s^{-l) s^{0) 1 
1 s^{0) sM 



J 



where 



sM = -X-V^{n)-VUn + l). 

Since K; assumes both and A on each cell [2k — 1, 2k], G Z, it must be the case 
that Su){2k — 1) = 0, i.e. s^) vanishes at all odd lattice sites. On the other hand, at 
even sites an inspection of cases shows that 

s^{2k) = \{uk+i - ojk)- 



B 



kk 
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Consequently, H^x is an infinite tridiagonal block matrix H^x = [Bjk]j,k& where 
each block Bji^ is a 2 x 2 matrix with Bji^ = I if |j — A;| = 1 and — 1 < j, A; < 1 and 

s^{2k) 
Su^{2k) 

Because of the 2x2 block structure, it is instructive to view H^^ x as a Jacobi-type 
operator on f{Z,C^)- That is, with u G fiZX"^), u{k) = {ui{k),U2{k)y , we write 

{H^,xu){k) = u{k - 1) + u{k + 1) + Bkku{k). 

Each Bkk niay be diagonalized via the transformation 

1 / 1 1 



W - 

which induces a unitary involution U on £^(Z, C^) given by 

{Uu){k) = Wu{k) for all A; gZ. 

H^^x is unitarily equivalent, via U, to another infinite tridiagonal block matrix M^^x = 
[Mjk]j^k& with each Mjk a 2 x 2 matrix; Mjk = I if |j — A;| = 1 and — 1 < j, < 1 
and Mkk = diag(s^(2A;), -s^{2k)). 

The fact that I and M^^ are diagonal means that M^^x decouples into a direct sum 
of "even" and "odd" parts. That is, if m G is expressed as the direct sum of 

its even and odd parts, corresponding to even and odd components, u = Uodd © ^even, 
then 

where are both discrete Schrodinger operators 

M^^ = -ho ±q^ = ho± q^, 
and the potential term is defined in terms of by 

qM = s^{2k) = X{uJk+i - Wfc). 

In light of the fact that 

the proof of Theorem 2^ reduces to showing 

ho±qu^> -V4 + A2 for all w G = {0,1}^. (34) 
If T : f2 — )■ f2 denotes the bijection defined by 

{Tu)k = -ujk + 1, (35) 

i.e. the 0-1-flip map, then 

qTcj = -q^ (36) 

Thus it suffices to establish (34) for ho + q^i. We do this by showing the existence of a 
positive function %p^, not necessarily belonging to ^^(Z), which satisfies the difference 
equation 

(^0 + g.)^. = -V4 + AV..- (37) 



This implies (34) by Theorem 3.3 
Let 



z±{X) ■= ^ 
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be the two distinct solutions of — z\/4: + + 1 = 0. Note that both are positive 
and z^{X)z^{X) = 1. We define ip^^ exphcitly by 

r z+{\)^U(^-^-'\ if fc>o, 

Mk) = < 1, if A; = 0, (38) 

[ z+{X)-^"^=>'+^^^'^^-^\ if A; < 0. 

Clearly, for all k E Z, tp^ik) > and 

Mk + 1) = Mk)z+{xf Mk - 1) = Mk)z+W-^''^''-'\ 

Thus 

-Mk - 1) + q^{k)i,^{k) - i)^{k + 1) 
= i^^ik) (-z+(A)2-'=+-i + \{uk+i - uju) - ^+(A)-(2-.-i)) 

= -v/4TA^^a.(fc), 

where the latter is easily verified separately for the four cases G {0, 1}, oJk+i G 



{0, 1}. This completes the proof of Theorem 2.6 
5.2. Density of States. 



□ 



Proof of Theorem 2^. Throughout the proof, we make use of the discrete Dirichlet 
Laplacian which on a set A C Z'^ is defined by 

where A^a is the edge counting operator defined in ([T]). 

If t < 0, then ^ is trivial. Thus, let t > 0. For any displacement configuration co, 
using the unitary involution (21) we have the unitary equivalence 

^A{L)(^) := ^oML) + ^^o,\a(l) = -ho^A(L) + -^K;|a(L) =: -H^L){uj). 

Thus 

#{keN:Ekih^(^L^{co))<E^{\)+t} = i^{keZ:Ek{Hj^^L)i^))>-E^{\)-t} 

= |A(L)|-n(i7^(^)(a;),-i5;_(A)-t) 
= |A(L)|-n«(^)(c^) + AlU(i), 

-^_(A) - t + A) 
= |A(L)| - n«(^)(a;) + AI|A(i)), E^X) - t) 
= |A(L)|-n«(^)(Ta;),E+(A)-t), (39) 

where, as before, in the last line = + XV. Also, n{A, E) denotes the number of 
eigenvalues of an operator A which are less than or equal to E and T is the bijection 
on Q defined by (35). 

If p = 1/2, then T is measure preserving on Q, 

F{T'^A) = F{A) (40) 

for all measurable sets A G Q, as it is induced by the measure preserving mapping 
a —a + 1 on {0, 1} with symmetric Bernoulli measure. 
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Dividing by |A(L)| and taking expectations in (39) gives 



^^E(|{A: G N : E.ihl^^^iu)) < E^X) + 1}\) = 1 - j^^E{g{Tuj)), (41) 

where g{u) = n{h^^^-^{uj), E+{X) — t). Since T is measure preserving, we have 
K{g{Tu)) = K{g{uj)). Now (41) immediately gives 



This gives (|6|) due to the independence of the IDS on the boundary condition. □ 



Proof of Theorem 2^. We give full details on how to do the proof at E^{X) and only 
give an outline of the modifications for proving the result at (7+ (A). 

Lower Bound at E-{X): Here we closely follow an argument developed in a similar 
context for the continuum random displacement model in [Ij and [3] . 

Let ^oA{L)('^) ^oA(L) + -^K;|a(l) for • G {D,N,P}. We will use the standard a 
priori bound 

NxiE) > ^^F^E.ihl^^^iu)) < E^ (42) 

see, for example, [121 (6.15)]. We will show that Ei{h^j^-^{u)) < E hj constructing a 
test function with Rayleigh quotient (\E'^, /i^.^s\E'^)/||\E'^|p < E. 



To this end, let 0i denote the strictly positive ground state of hjl^ + X62 normalized 
so that 0i(l) = 1, and let denote the strictly positive ground state of h^_^ + X6i 
normalized so that 0-i(l) = 1. Using uniqueness of the positive ground state, we 
have 0_i(2) = 0i(2)~^ In fact, (pti given by 4>Li{n) = 0i(2 — n + 1) is a positive 
ground state of + X6i. For ease of notation, we denote r := 0i(2). We know by 



Lemma 4^ that r 7^ 1 and will now assume that r > 1. IfO<r<l, we can do the 
following construction from "right to left", choosing the test function "^^^ such that 
*^(2L) = 1. 

Given u, we construct "if^ by concatenating rescaled versions of 0_i and 0i. The 
test function "^i^ is defined piecewise on A^, 1 < < L, by 

^.(n) = r^t>'("Va,H(^ - 2{k - 1)) n G A^, (43) 

where ak{uj) = — 1 if = and ak{uj) = 1 if = 1. This choice of rescaling 
guarantees that '^^{2k) = \E'^(2A; + 1) for all k and thus, by the properties of Neumann 
boundary conditions on all A^, that h^^-^{ijj)'^i^ = _Eo(1)^I'lj- However, as we have 



shown in the proof of Theorem 2.2, Eo{l) = E^ = mma{ho + K;*). Thus we have 



^oA(L)'^'^ = -E-^&o; and we may write 



The Dirichlet and Neumann operators only differ at the endpoints of A(L), so most 



terms in the numerator of (44) cancel; what rests is 



(^^.^A(L)M^^) _ 2(l + v^^(2L)^) 
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With (42), we have 




2(1 + ^^(2L) 



2^ 



<E-E_ 



It follows from the definition of "^^^ that 

log^<.(2j)/^.(2j - 1) log^^(2j + l)/^.(2j - i: 



log r log r 

(where we use that '^^{2j) = \l'a;(2j + 1)) and 

^^(2j)2 = e'^^^°s^ (46) 

where Sj := ai + 02 -\ haj. As F{ak{uj) = 1) = f{uJk = 1) = 1/2 and P(afc(w) = 

—1) = F{uJk = 0) = 1/2, the aj{u) are independent symmetric Bernoulli random 
variables with values in {±1}. Therefore, the process Sj, j = 1,2, . . ., is a simple, 
symmetric random walk. If Y := maxj=i ^Sj, then it is a consequence of the 
reflection principle for symmetric random walks ([8J) that 

P(r >Vl\Sl<0) = P(5l > 2/L). (47) 

The latter converges to 

POO 



as L — !■ 00 by the central limit theorem. 

Denote by Al the event {uj\Y > Vl and Sl < 0}. If F > /L, then ^^^^ ^Mjf > 
g2A/Liogr_ rjnj_^g condition 5*^ < means that equal or more single site potentials sit at 
the left than sit at the right on A(L). As r > 1, it is clear from (46) that \E'a;(2L)^ < 1 
and we have 



p(^Ei(/.^(^)(c^))<i?^ > P 



2(1 + ^^(2L)- 



<E-E\A, 



> F(^4e-^^^°^'' < E - E_\Al^F{Al 



P(Ai) 

F{Y >\^\Sl< 0) ■ F{Sl < 0) > > (49) 



if 4e-2v^i°g^ < E - E^ and L > Lo (where we have used (|47| and (|48|)). If E is so 
close to E_ that 



22 ROGER NICHOLS AND GUNTER STOLZ 

then we can choose L > Lq for which 

4g-2Vi:3Tiogr >E-E^> 4e-2^i°g'-. 

Therefore there are constants Ci, C2 > such that 



CiL < 



logr 



<C2L. 



(50) 



(51) 



From (49) we get A^a(^) > ^7, so that 



Nx{E) > 



logr 



Cn 



2C,\\ogl{E-E^: 



> 



C 



\og\E-E^] 



for E — E^ sufficiently small. 

Discussion of Modifications for Lower Bound at G+(A): We want to consider a 
lower bound for N\{Gj^{\) + e), with e > 0. The idea is to express the integrated 
density of states A^a at the gap edge 6*+ (A) in terms of th e integrated density of states 
of the operators ± from the proof of Theorem 2.6, call them N±, at — a/4 + A^. 
One then estimates N± from below. 

Let L G N and define 



- 2 + 



A^ 



(52) 



It follows from Theorems 

tL 



2.3 



and 



2.6 



that hf[ i^) eigenvalues in the interval 



{G-{X),G-i-{X)). Thus, H^^^ > — v4 + A^. Moreover, via the same construction as in 
Section 



5.1 



^^ is unitarily equivalent to the direct sum 



Ht,X = (^o![l,2L] - Qi^\[l,2L]) © {^0,[1,2L] + ?a;|[l,2L]) =: Jl 

Therefore, setting r]{e) := 2e(G+(A) - A/2) + e^. 



e - 



Nx{G+{X) + e)-Nx{G_{X)-e) 



L-5- 



jim ^E(#{fc : i?.(/iJ,2L]M) e [G_(A) - e, G+(A) + e]} 

lim ^e(^#{ A; : i?,(J^ © J+) < -v/4TA^ + r7(e)}^ 
hm ^Ef#{A; : E,(J^) < -./ITX^ + r]{e)} + 
: Efc(J+) < -v/4TP + r/(e)}^ 



-iV_(-y4TA2 + r7(e)). 
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The last equality uses that N_ = N^, which follows from the fact that T defined by 



(35) is measure-preserving on Q and satisfies qxai = —Qui- In particular, 



iV,(G+(A) + - ^ > (-v/4TA^ + 26(G+(A) - . (53) 

One can prove that has a 1/ log^-singularity near + with essentially the 
same techniques used for Nx near E^{X). For one chooses a different trial function 
\E'a;- The appropriate choice is "^uj '■= V^wlfi.L], where ip^j is the positive function defined 
by (38) in the proof of Theorem 2.6 With this choice, is, up to an error term, an 
eigenfunction of /i^yv(L) ~ Q^^Iml) with eigenvalue —y/A + A^. More precisely, 

(^a.,K[l,L]-g.|[l,L])*.) 

= -VI+X^ ll^^f + ^^(1)2 - ^^(1)^^(2) - ^.(L)^^(L - 1) + 
Thus we have the lower bound 

[ly - - ^.im.jL - 1) + ^ ^ ^ v/4Ta^^ 

From here on the proof is completed in a similar fashion to what was done above 




for E^{X). Key is the simple multiplicative structure of ip^i as defined in (38), which 



leads to considering the symmetric random walk defined by the Bernoulli variables 

log^^(A: + l)/^^(A;) 

ttfc := ] 7TT = 2WA:+i - 1. 

log Z+{X) 

If A > 0, then 2+ (A) > 1 and the proof goes through with the same argument as 
above. If A < 0, then < 2;+ (A) < 1 and one can work from "right to left", similar to 
what was indicated for the case < r < 1 above. 

□ 

6. Concluding Remarks 

(i) Our results for d > 1 fall short of what was proven in [2] and [3] for the continuum 
displacement model in several respects: 



In part (i) of Theorem 2.1 we have required sign-definiteness of q. The corre- 
sponding result for the continuum from [2] only requires that Eo{a) does not 
vanish identically in a. In other words, in the continuum the multi-dimensional 



analogue of part (ii) of Theorem 2.1 holds (noting that is the spectral mini- 
mum of the continuum Neumann Laplacian, while —2 is the spectral minimum 
of the discrete one-dimensional Neumann Laplacian). 

We have yet to prove any uniqueness results for the set of periodic configu- 
rations u with min(T(/i^) = Emin in the case of d>2, but we conjecture the 
following analogue of a result which was shown to hold for the continuum case 
in|3|: 



24 



ROGER NICHOLS AND GUNTER STOLZ 



Conjecture 6.1. If d > 2 and q ^ is sign- definite, then u* is (up to 
translation) the unique periodic configuration with min cr(/;,(^.) = Emm O'nd 
maxcr(/itj.) = E„ 



-'max ■ 



Both of these shortcomings of our results are due to the lack of unique continuation 
properties of the discrete Schrodinger equation, which were used in this context in 
the continuum in [21 E] . 

(ii) It has recently been shown in [15] that the continuum multi-dimensional random 
displacement model, under symmetry assumptions similar to the ones used here, ex- 
hibits Anderson localization at energies near the bottom of the almost sure spectrum. 
Proving this also in the discrete case remains a challenging problem. A uniqueness 
result as conjectured in (i) above would allow to prove Lifshitz tails for the IDS with 
a strategy developed in [T7| (and used in [IS]). However, due to the "discrete" na- 
ture of the randomness we can not expect a Wegner estimate to hold for the discrete 
displacement model. The new type of multiscale analysis recently developed in |5] to 
circumvent this problem for the continuum Bernoulli- Anderson model does not work 
on the lattice, again due to the lack of unique continuation properties. 

(iii) For the one-dimensional Bernoulli displacement model we have no result like 



Corollary 2.7 if |A| > 2, but our conjecture is 
Conjecture 6.2. For any A 7^ 0, 

SA = ^T(/l^,,A)Ua(/l^i,A), (54) 

where is the displacement configuration with components uj\ = 1 for all k E Tj. 

If Conjecture[2]is true, it would mean that consists of exactly six bands separated 
by five non-vanishing spectral gaps. Other than the fact that this would be a natural 



generalization of what we can prove for |A| < 2 (where Corollary 2.7 followed as 



a consequence of (54)), we have numerical evidence for this conjecture which was 
provided to us by Jeff Baker. Figure |2] shows numerical plots of the density of states 
(through histograms of the eigenvalues of finite volume restrictions averaged over 
multiple sample configurations) as well as the IDS for A = 3. The same pattern arises 
numerically for other values of A > 2. 



iv) Note that Figure M also illustrates the 1 / log -lower bounds for the growth of 



the IDS found at the band edges E±{X) and G±{X) in Theorem 2.9 One of the reasons 
which make this bound interesting is that it was proven in |6j that for general ergodic 
Jacobi matrices the IDS satisfies an upper bound of the form \N{E) — N{E')\ < 
C/\log\E — E'W locally at all energies, i.e. that N{E) is log-Holder continuous. It 
would be interesting to find out if our lower bound is sharp, i.e. if there is also an 
upper bound of type 1/ log^. 

Known examples showing the optimality of the 1/| log | upper bound were provided 
in [6] and, more recently, in [9]. These examples are in terms of quasi-periodic and 
limit-periodic potentials. As opposed to our example, these are non-random in the 
sense of having long-range spatial correlations. 

(v) The eight additional band edges which would exist for |A| > 2 if Conjecture |2] 
were to hold show different numerical properties in Figure [2j This is most drastically 
apparent in the density of states plot, where delta-peaks appear at the edges addressed 



by Theorem 2.9, but not at the other eight edges. Based on the numerics we state 
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Density of States 




Integrated Density of States 




Figure 2. The density of states and integrated density of states for 
the one-dimensional BDM with A = 3 and p = 1/2. 



Conjecture 6.3. // |A| > 2, then the IDS has thin tails (e.g. of Lifshits type) at the 
eight other band edges ofIl\. 
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